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In this note we report on geometric scaling in inclusive eA scattering data from the NMC and 
E665 experiments. We show that this scaling, as well as nuclear shadowing, is expected in the 
framework of nonlinear pQCD at small x based on a simple rescaling argument for ep scattering. 



INTRODUCTION 

Perturbative QCD, in particular where based on fac- 
torization theorems for high energy processes, has been 
very successful in describing data from high energy ex- 
periments. However, perturbative evolution has almost 
always been linear in the parton densities, which was, 
and very often still is, quite sufficient to describe high 
energy scattering. 

There are, however, situations where this is bound to 
fail. Among those are high energy scatterings at energies 
large enough to create a gluon medium in which multi 
gluon correlations can no longer be neglected. This has 
first been discussed in the context of deep inelastic scat- 
tering of electrons off (large) nuclei, where both the large 
number of nucleons and the high energy lead to an in- 
crease in the number of gluons involved in the scatter- 
ing |, |, 1 1110 I, §• In this regime nonlinear 
effects, leading to a taming or saturation/unitarization 
of F2 , should be included. 

In the guise of a very successful phenomenological fit, 
this idea has even been applied to the ep experiment at 
HERA with its strong and, if untamed, unitarity violat- 
ing rise of the proton structure function F% (xbj, Q 2 ) with 
decreasing [0 • This has led to the observation of 
"geometric scaling" satisfied by HERA data in the small 
Xbj region Q. Our main focus here will be an attempt 
to extend this geometric scaling idea to nuclei and to 
check if it is already confirmed by available eA data. 

The question of the onset of this satura- 
tion/unitarization in inclusive QCD observables with 
increasing energy (decreasing Xbj), in particular in the 
proton structure function (sbj, Q 2 ), has been the 
subject of active discussions for many years (see a de- 
tailed discussion of this subject in H, |l2| and references 
therein.). Phenomenologically, the key feature at small 
Xbj is the simple structure of the total j*p cross section 
or F^' A that arises directly from the underlying physics 
picture. Viewed in the infinite momentum frame of the 
target, the virtual photon splits into a qq color-dipole 
(treated as eikonalized Wilson-lines) with a distribution 
characterized by Q 2 . This dipole then punches through, 
and interacts with the target, which it sees as a pancake 
of infinitesimal longitudinal thickness. Consequently the 



cross section appears as a convolution of the square of 
a photon wave function, which gives the probability to 
create the qq pair, and a dipole cross section, which 
contains all the information about the target and the 
strong interaction physics. 

The latter carries x^j dependence which arises from 
gluonic fluctuations which induce an increase in the num- 
ber of gluons ]l9| ]. From this point of view, the target 
will display a growing density of gluons per transverse 
area (integrated over the longitudinal extent of the tar- 
get) that drives the system into a saturation regime. The 
corresponding formula for the total cross section in "f*p 
is then 
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where a is the longitudinal momentum fraction of the 
q or q, and z their relative transverse separation. Note 
that the Q 2 and ibj-dependence is clearly separated into 
wave function and dipole cross section. This reflects 
the natural distinction between the transverse and lon- 
gitudinal directions seen at high energies. Q 2 sets the 
scale for transverse resolution, but longitudinal informa- 
tion about the target is subsumed in the Xbj dependence 
of the dipole cross section. The idea is that the dipole 
cross section depends on x^j only via a scale Q s {xbj) as 
CTdi P oicObj,z 2 ) = o-dipoie^QsO^bj))- Q s is the satura- 
tion scale, i. e. the scale below which the gluonic content 
of the target becomes "black" or Odipoie reaches its large 
size asymptotics. 

Together with Eq. (|l|) this implies that <7tot and hence 
F% scales according to 

o-tot(^bj,<3 2 ) = ertot(xo, Q2( Xh .} Qo) ■ ( 2 ) 

This phenomenon is called geometric scaling and, as 
already mentioned, was first observed to be present 
in HERA data @. The fit used in |t2| assumes 
an "eikonal" shape for erdipoic that saturates for large 
dipoles. As Xbj decreases this saturation affects smaller 
and smaller distances, inducing Q s to grow, a fact which 
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was parametrized by 



.(Xbj) 
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The theoretical basis for such saturation/unitarization 
features, with a scaling behavior as in Eq. (Eh, together 
with precisely the functional dependence of Q s as in 
Eq. ([|), lies in nonlinear evolution equations that try 
to predict this Xbj dependence from first principles. 

The most promising approach to properly include sat- 
uration/unitarization is the JIMWLK equation which 
resums the leading ln(l/iEbj) terms in all N-point cor- 
relation functions of the participating fields not just the 
leading correlators as in fc_L-factorization. This leads to a 
non-linear renormalization group equation (RGE) for the 
generating functional for these N-point correlators which 
slows down the rapid growth of structure functions at 
small Xbj while at the same time preventing the system 
from drifting into the infrared. This is the basis for a 
theoretically sound, selfconsistent treatment. Lineariz- 
ing the JIMWLK equation, one recovers the well-known 
BFKL equation (see |7|, [l| for details) as its low den- 
sity limit. We thus learn to interpret the loss of infrared 
safety in BFKL as being due to the absence of any non- 
linear corrections and the predicted too rapid a growth 
in F^XbjjQ 2 ) as signals of BFKL having exceeded its 
range of validity. 

The full content of the JIMWLK equation can be 
accessed numerically: being a functional equation of 
Fokker-Plank type it can be rewritten as a Langevin 
equation with white noise. This then can be discretized 
and solved using lattice methods juj. This yields the 
Xhj dependence of Cdipoie and in particular the evolution 
"rate" A of Eq. ©. 

Of particular interest to us in this context is the fact 
that the evolution equation is target independent. Tar- 
get dependence arises only from the initial conditions for 
the Xbj dependence of the dipole cross section (related to 
the initial distribution of gluons in the transverse plane), 
which will, in particular, be A dependent. These facts 
will allow us to compare ep and eA scattering and in- 
vestigate whether one of the basic predictions of this ap- 
proach, geometric scaling, can also be found in inclusive 
eA processes. 



GEOMETRIC SCALING AND NUCLEAR 
SHADOWING WITHIN NONLINEAR PQCD. 

Before talking about geometric scaling in eA scatter- 
ing, one should address nuclear shadowing or the suppres- 
sion of \f£ to F% at small Xbj- If this is at all possible 
it increases our phenomenological leverage considerably 
by allowing us to compare experiments with different A 
and Xbj directly. This is of particular importance as the 



available eA data do not offer any leverarm in Xbj for a 
given Q 2 or vice versa. 

In fact, the leading effect of nuclear shadowing can 
be easily incorporated in this approach by realizing that 
the change in transverse size simply changes the overall 
normalization of <7di po ie while the change in longitudi- 
nal extent can be incorporated by increasing our initial 
Q s by a factor A s . In the simplest case, assuming the 
distribution of spectator partons in the target to be ho- 
mogeneous, S = 1/3. Eq. for F2 ~ a ■ Q 2 then turns 
into 



x bi \ 2X F 2 A (x hj ,Q 2 

x J ,42/3+1/3 
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Note that < Xbj < A since Xbj is given per nucleon. 
The power 2/3 in the overall A dependence stems from 
a simple geometric factor, the area of the target, in the 
impact parameter integral in Odipoie- Thus we have an 
unambiguous prediction for -jF^ which states that not 
only can the observed nuclear shadowing in the eA data 
be explained by a simple rescaling of the variable Q 2 
with A 1 / 3 , but also that all data should, plotted vs. 

T = (~c?f) 3*73' ^ e on ^ e same curve as the data for 
Ff • Hence, this approach also predicts geometric scaling 
in inclusive eA scattering analogous to the one observed 
in ep scattering. 

GEOMETRIC SCALING IN \f£ 

Note that Eq. (0) is strictly true only in the asymptotic 
region of small Xbj, large Q 2 and large nuclei. Therefore, 
we expect that there might be nonperturbative correc- 
tions which can spoil the form of Eq. (0) and consequently 
we perform a fit to the combined nuclear data with the 
following form: 
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where 7 and 5 are the fit parameters, which ideally would 
be and 1/3 respectively. 

The starting point of our analysis is the published 
NMC Q and E665 ||, data for the ratios of nuclear 
to deuterium structure functions 2F/ / ' AF® . In Fig. |l| 
we show the ratios as functions of x. From these ratios 
we obtain the nuclear structure functions F 2 as follows: 
we first note that NMC data show the deuterium to pro- 
ton structure function ratio being F® /2F% — 1 to better 
than 95% accuracy in the dynamical range of interest 
(xbj < 0.1) which is well within the error margins of 
the nuclear F2 measurements. Thus, we convert F^/F^ 
to structure functions F^ by multiplying the ratios by 
F%, using the geometric scaling ansatz to extrapolate F% 
to the appropriate values of Xbj and Q 2 . 
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FIG. 1: Ratio of 2F^/AF? vs. t for NMC and E665. Note 
the increasing steepness with A within the NMC and E665 
datasets respectively. 
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FIG. 2: Comparative phase space of HERA, NMC and E665. 



This rescaling is necessary since the kinematic ranges 
of the HERA, NMC and E665 experiments are not iden- 
tical as shown in Fig. |[ 

Now we are ready to use Eq. (JsJ) to compare nuclear 
and proton structure functions. In Fig. [I] the rescaled F£~ 
is seen to fall on top of the dashed line representing the 
geometric scaling fit to the HERA F$. For comparison we 
also show the quality of the scaling fit within the HERA 



FIG. 3: Scaling behavior of NMC and E665 F 2 A data vs. r = 

. \2A 2 

— I .i i, ■ The vertical axis corresponds to the l.h.s. of 

x ° L A ' 

Eq. (0) . The dashed line corresponds to the geometric scaling 
curve obtained from HERA data. These are shown offset by 
a factor of 5. 



data, offset by a factor of 5. 

The best fit to the combined nuclear data is achieved 
for 7 = 0.09 and 5=1/4 which are close to but not quite 
the asymptotic values the nonlinear pQCD approach pre- 
dicts. The remaining parameter A ~ 0.18 in Eq. (|^) 
is already determined from the geometric scaling fit to 
the HERA data. With this value for A, the saturation 
scale at the upper end of the HERA shadowing region at 
x = 10 _1 comes out to be marginally "perturbative" , 
Ql = 1 GeV 2 . 

We found that the NMC calcium data allows a smaller 
value of 7 and, more importantly, a value of S = 1/3 if 
the data on lighter nuclei is left out. This can be seen 
best in a log-linear plot of the scaled F2 ratios shown in 
Fig. [I], where the plateau from r = 0.1 to r = 10 shows 
good overall scaling behavior. 6 = 1/3 would lift the 
NMC Ca data fully up to 1. This seems to indicate that 
the Ca nucleus is large enough for the basic assumptions 
of the nonlinear pQCD approach to be valid. A very 
encouraging sign, indeed. 

Also note that there is a consistent difference in the 
E665 and NMC data with the E665 data being much 
flatter in the ratio 2F^ /AF® , as evidenced in the com- 
parison between the two calcium data sets and shown in 
Fig. |1| This apparent inconsistency between the two data 
sets disappears when the ratio YIF^ /AF^ is compared 
between the two experiments |i"8] . This question needs 
to be addressed experimentally. 
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FIG. 4: Plot of scaled NMC and E665 F^/AF% data vs. 

T = (^f) 41/3 • This corresponds to "dividing the data 
shown in Fig. pa by the dashed line." 



Notwithstanding all of the above, nuclear data from 
the proposed Electron Ion Collider at BNL on heavy nu- 
clei at small and large Q 2 is necessary to be absolutely 
definitive about geometric scaling in inclusive eA scat- 
tering. These experiments will by design provide data at 
smaller x\,j but one should take care to open up the phase 
space region. Only if we have access to a range of Xbj 
values at any given Q 2 (and vice versa) for each species 
A, can we obtain independent fits for the evolution rate 
A. This would be a prerequisite to experimentally dis- 
entangle and confirm the A and x dependence extracted 
simultaneously above. As the errors on the nuclear data 
get smaller one would also require a direct measurement 
of F^/AF^ that would allow us to do away with our 
rather cavalier treatment of F 2 I) /2i 7 '|\ 



CONCLUSIONS 

To summarize, we demonstrate for the first time that 
the NMC and E665 data for the nuclear structure func- 
tion -jF^ also exhibit, modulo the caveats given above, 
geometric scaling as already seen in the nucleon structure 
function F%. This has been confirmed using the shadow- 
ing features of the nonlinear pQCD approach, which, as 



the dominant effect, consists of a simple rescaling of Q 2 
by A 1 / 3 . With the presently available nuclear data which 
are at rather low Q 2 we had to allow slight modifications 
in the powers entering the the A dependence in order to 
get a perfect fit. The NMC data for calcium, the heav- 
iest nucleus available within this dataset, would in fact 
be compatible with the A 1 / 3 scaling. Given the limited 
phase space regions covered in the data sets we find it 
very encouraging how well the geometric scaling in eA 
follows the one in ep. 
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DFG, AS by BMBF and HW by the DFG Habilitanden- 
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